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Using different techniques, and Fermi-liquid relationships, we calculate the variation with applied 
magnetic field (up to second order) of the zero-temperature equilibrium conductance through a 
quantum dot described by the impurity Anderson model. We focus on the strong-coupling limit 
U ^ A where U is the Coulomb repulsion and A is half the resonant-level width, and consider 
several values of the dot level energy Ed, ranging from the Kondo regime ep — Ed ^ A to the 
intermediate-valence regime ep — Ed ^ A, where ep is the Fermi energy. We have mainly used 
density-matrix renormalization group (DMRG) and numerical renormalization group (NRG) com¬ 
bined with renormalized perturbation theory (RPT). Results for the dot occupancy and magnetic 
susceptibility from DMRG and NRG+RPT are compared with the corresponding Bethe ansatz re¬ 
sults for U —>■ oo, showing an excellent agreement once Ed is renormalized by a constant Haldane 
shift. For U < 3A a simple perturbative approach in U agrees very well with the other methods. 

The conductance decreases with applied magnetic field for dot occupancies ~ 1 and increases 
for Ud ~ 0.5 or Ud ~ 1.5 regardless of the value of U. We also relate the energy scale for the 
magnetic-field dependence of the conductance with the width of low energy peak in the spectral 
density of the dot. 

PACS numbers: 75.20.Hr, 71.27.-l-a, 72.15.Qm, 73.63.Kv 


I. INTRODUCTION 

In the last years an enormous amount of research in the 
field of nanoscience has been devoted to the transport 
through semiconductingir— and moleculari^rJJ quantum 
dots (QDs), and to manifestations of the Kondo effect 
in these systems. The semiconducting QDs are artifi¬ 
cial atoms created in two-dimensional electron gases by 
a suitable application of electrostatic voltages and are 
characterized by a high tunability of the parameters. By 
contrast in the molecular QDs, the molecule itself can be 
changed. 

Many of these systems are described by the Ander¬ 
son model for a magnetic impurity with spin 1/2 [Eq. [T] 
below], with constant hybridization V between the im¬ 
purity level and the conduction electrons, whose unper¬ 
turbed density of states p, can also be assumed constant. 
The main parameters of this model can be taken as the 
energy Ed of the level localized in the QD relative to the 
Fermi energy, which we take as = 0, half the resonant 
level width A = npV^ (in the literature sometimes the 
total width T = 2A is used), and the Coulomb repulsion 
U. As we shall see, the conduction electron band width 
2D also plays a role, although a minor one. In the Kondo 
regime, —Ed A and Ed+U A, all quantities depend 
on a single energy scale Tk ^ D^/pJex^p[—l/{pJ)] with 
J = 2V‘^U /[—Ed[Ed-\-U)\^ and the localized spin at the 
QD is compensated by the conduction electrons result¬ 


ing in a singlet ground state and maximum conductance 
(unitary limit) at temperatures T <C Tr-.— 

Recent experiments have studied the scaling laws for 
the conductance through one QD in the Kondo regime 
for small (<C Tk) bias voltage 14, temperature T and 
applied magnetic field This stimulated further the¬ 
oretical work on the subject concentrated mainly on the 
effect of a finite (although small) bias voltage, which is 
a tough non-equilibrium problemii2ri22 While it would 
be desirable to express the non-equilibrium properties 
in terms of thermodynamic ones, this task seems pos¬ 
sible only for a limited number of coefficients in the ex¬ 
pansion of the conductance!^ Instead, at T = 0 and 
at equilibrium (14 = 0), the magnetic-field dependence 
of the conductance, characterized by a coefficient cb 
which was addressed theoretically recently;^^— can be 
expressed in terms of the magnetic susceptibility and the 
second derivative of the dot occupancy with respect to 
the magnetic fieldr^ using an extension of the Friedel 
sum rule^ to finite magnetic field for a spin conserv¬ 
ing impurity model ;^^'^^ used before in magnetotransport 
through quantum dotsi^^— 

The effects of magnetic field in systems of one and sev¬ 
eral QDs have been studied before experimentalljAi^Zr— 
and theoretically2^i2^i^“— using several methods. The 
coefficient cb has been calculated using numerical renor¬ 
malization group (NRG) and a “superperturbation” 
approach2^i21 recently i^^i^° Unfortunately, the results 
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presented in Ref. l28l . while exact at the symmetric 
point, are in general inaccurate because they miss a 
term proportional to the second derivative of the occu¬ 
pancy with respect to the magnetic field^S (the effect of 
this term is responsible for the difference between top 
and middle lines in Fig. [T] below ). Correct results 
were later presented (22 but only in the weak coupling 
regime U < 3A. For U —>■ oo, the coefficient cb was 
calculated using a slave-boson mean-field approximation 
(SBMFA)/— However, direct comparison with NRG re¬ 
sults shows that the SBMFA does not provide a reliable 
dot occupancy nd^ and since Ud and its second deriva¬ 
tive with respect to the magnetic field enter the expres¬ 
sion of cs [see Eq. dSD], the results of Ref. [2^ are only 
qualitatively valid. 

On the other hand, for U < 3A, the possibility to 
reach the Kondo regime {—Ed ^ A and Ed + U A) 
is questionable and limited in any case to the symmetric 
situation Ed = —U/2. In general, in the Kondo regime, 
the spectral density presents two charge-transfer peaks 
at Ed and Ed + U of total width 4A and a Kondo peak 
at the Fermi energy of total width ^ 2Tk^— In the 
symmetric case Ed = —U12 only one peak (at the Fermi 
level) is present for U < 3A)^ and as we shall discuss in 
this work, its half width at half maximum for U = 3A is 
0.54A. The similarity of all energy scales and in particu¬ 
lar the large Tk implies that one has to reach magnetic- 
field energies gfisB of the order of U to obtain a splitting 
in the spectral density twice which is a distinc¬ 
tive signature of Kondo physics For smaller mag¬ 

netic fields the splitting is smaller. A splitting roughly 
consistent with 2gg,BB was observed in nonequilibrium 
transport through QDsiii^Zr— More recent experiments 
obtain a somewhat larger splitting. This is probably due 
to the fact that these experiments should not be inter¬ 
preted in terms of the equilibrium spectral density and 
non-equilibrium calculations are necessaryi^^f^^— In any 
case, to detect clearly the splitting one needs devices for 
which Tk U. 

Experimentally in semiconductor QDs a wide range 
of values of [//A are possible. Typical values of U 
are around 1-2 meVfii^ii and A can vary between 10 
to 200 p.eV, leading to ratios [//A near 10 in the first 
experiments^^ 2 and 4.5 in the two devices studied by 
Kretinin et al^, and larger than 50 in a recent system of 
two QDs^ In contrast in molecular QDs the expected or¬ 
der of magnitude of [/ is ~ 1 eV for small molecule o^^d^ 
or ^ 0.1 eV for large molecules 3, while typical values of 
A are of the order of ^ 1 meVj ^^d^ Thus, the ratio [//A 
is several orders of magnitude and in practice one can 
take [/ —>■ oo. These facts indicate the need to extend 
the reliable calculations of Cs (so far limited to the weak 
coupling regime U < 3A) to more realistic values. 

In addition, the molecular QDs are characterized by a 
high asymmetry of the coupling of the QD to the leads, 
and as a consequence, the shape of the differential con¬ 
ductance G = dl/dVb near zero applied bias voltage 14 
reproduces the spectral density of the Kondo peakj ^"^4° 


This renders possible to relate the energy scale of cb with 
the width of the measured zero-bias anomaly as we shall 
discuss. 

In this work we calculate cb as a function of Ed for sev¬ 
eral values of U and in particular in the strong-coupling 
limit C/ —>■ oo. We also calculate the half width at half 
maximum Ap of the low-energy peak in the spectral den¬ 
sity of the dot for for several values of Ed and U. As Ed 
increases from the symmetric case Ed = —U/2, cb de¬ 
creases, changes sign in the intermediate-valence regime 
Ed ^ 0 (where the occupancy Ud 0.6) and becomes 
large and negative in the empty-orbital regime Ed > 0. 
As a function of n^, cb looks qualitatively similar for all 
values of U. 

The paper is organized as follows. In Section II we 
present the model and equations that relate cb with ther¬ 
modynamic quantities. In Section III we describe briefly 
the different methods used. Section IV contains the re¬ 
sults and Section V is a summary. 

II. MODEL AND FORMALISM 

The model is an Anderson impurity one in which a 
single localized level in a QD is hybridized with two con¬ 
ducting leads, 

H = X/i/feo- + X/cr BdTLda + Und\TLd\, 

+ 'Lukc^iYi'kdlcka- + H.C.) - gUBBs^, (1) 

where cj)^^creates a conduction electron at the lead v 
with momentum k and spin cr in the conduction band, 
dj. creates an electron in a localized level of the QD, 
ndcr = d\da, and = {nd\ - nd^)/2. By consider¬ 
ing appropriate linear combinations of the electrons of 
both leads, the model is mapped into a single-channel 
model with resonant-level half-with at half-maximum 
A = ~ Efe)? which we assume indepen¬ 

dent of energy w. 

At zero temperature (T = 0) and at equilibrium (bias 
voltage 14 = 0), the contribution to the conductance of 
each spin for a given magnetic field Ga{B) is propor¬ 
tional to the corresponding density of states of the dot 
level pda{i^,B) at the Eermi energy w = Oi^i^ In turn, 
this quantity is related to the occupancy for the corre¬ 
sponding spin by the Friedel sum rule^ generalized to 
finite 

/n I:.^ sin^ (imda) 

Pda{0, B) = - — -, (2) 

where for simplicity we denote as nda the expectation 
value of the corresponding operator. This allows to ex¬ 
press the change of conductance by an applied magnetic 
field in terms of the occupancies 

G{B) ^ pd^{f),B)+pdx{f),B) 

G(0) Pdt(0,0)+pd^(0,0) • 







3 


Expanding Uda up to second order in B one 
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XB d'^Ud B'^ 
dB ^ 4 


+ 0(i?3), 


( 4 ) 


where = nd-\ + nd\^, X is the magnetic susceptibility, 
a = 1 (-1) for spin up (down) and the quantities in the 
second member except B are evaluated at B = 0. 

From these equations, and defining cb and Tq (an en¬ 
ergy scale of the order of Tk) by2^“— 


G{B) _ _ (ggBBV 
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III. METHODS 

As explained in the previous section, the calculation of 
the coefficient cb that describes the magnetic-field depen¬ 
dence of the conductance at equilibrium and zero tem¬ 
perature, reduces to the calculation of the dot occupancy 
Ud, the magnetic susceptibility x and the second deriva¬ 
tive of the occupancy with magnetic field d^Ud/dB^. We 
have used four methods to calculate these quantities. 


A. Density-matrix renormalization group 

We used density-matrix renormalization group 
(DMRG)^ to solve the impurity Anderson model for 
B/A = 3, 8 and also for the infinite U case. The 
band was discretized using a Wilson chaliJ^, with a 
discretization parameter A. This presents benefits that 
have already been introduced in the literature<^“— In 
the range of parameters considered, the results do not 
change significantly when increasing the number of sites 
N beyond a certain value A*(A), as we have numerically 
verified for several values of A. We found that setting 
A = 4, A = 50, and retaining m = 600 states was 
enough to assure convergence, being the truncation error 
(the weight of the neglected states in the density matrix) 
10“® in the worst case, reassuring the reliability of the 
calculation. The DMRG results in this paper correspond 
to these values of A, A and m. 

As a consequence of this discretization, the density of 
conduction states at the Fermi level pa(0) decreases with 


respect to the continuum limit A —> 1. We have calcu¬ 
lated pa(0) numerically for the Wilson chain without dot, 
and determined the hybridization Va of the dot with the 
first site of the chain from the condition 7rpA(0)V^ = A. 

We have chosen a half band width D — lOOA since the 
ratio D/A is very large in real systems. 

We have calculated the ground state energy and the 
occupancy of the impurity for different values of the mag¬ 
netic field B, which was applied to the impurity site only, 
as in Eq. O- For the calculation of the susceptibility 
and d'^Ud/dB^, the general criterion used was to calcu¬ 
late the ground-state energy and magnetization at the 
dot for ten different values of the magnetic field B such 
that Bjjiiji ^ B <i Bfjidx-i where gpBBmin = 0.005T^, 
and Bmax ~ ISOBmm- The curves were smooth and 
very well approximated by linear or quadratic polynomi¬ 
als, depending on the case. The susceptibility was calcu¬ 
lated both as the derivative of the magnetization of the 
dot with respect to B fitting the curve with a straight 
line through the origin, and as the second derivative of 
the ground state energy with respect to B fitting with a 
quartic polynomial and extracting the coefficient of the 
quadratic term. In all cases the agreement was complete. 


B. Numerical renormalization group plus 
renormalized perturbation theory 

We have used the standard Numerical renormaliza¬ 
tion group (NRG)^iS to calculate the occupancy Ud 
and its second derivative with respect to the magnetic 
field d’^UdjdB"^. The latter was obtained calculating 
Ud for nearly 10 different magnetic fields in the inter¬ 
val 0 < gpBB < A, where A is defined below and fitting 
the points with a B^ dependence. 

We used a discretization parameter A = 3.5 and trun¬ 
cated the spectrum, after the fifth iteration, keeping 
up to 2000 states. In contrast to the DMRG calcula¬ 
tions, we take D = lOA and do the numerical calcula¬ 
tions using a hybridization Va = y/A\V between the dot 
and the first conduction site in the Wilson chain, where 
A = ttV'^/{2D) and 


Aa 


A + 1 
2A-2 


InA. 
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This expression has been suggested to obtain the correct 
Kondo temperature in the Kondo limit i^^i^^ . 

For the magnetic susceptibility y one has the problem 
that in the Kondo regime (large U and A <C \Ed\,Ed-\-U) 
it oscillates as the iterations increase (lowering the tem¬ 
perature). One way to solve this problem is to use A ^ 1 
and average over several realizations of the logarithmic 
grid)^ but we obtained only a moderate improvement. A 
better method is to use the full density matrix within the 
NRG, although still in the strong coupling limit, values 
of the Wilson ratio R above 2 and about 3% larger than 
the exact ones were obtainedi^ Here we have obtained 
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renormalized parameters A and U that describe the 
low-energy physics, following the procedure explained by 
Hewson et al. in Ref. llS The only difference is that we 
interpret that in Eq. (42) of that paper, the first mem¬ 
ber refers to A a the renormalized A for A 7 ^ 1, which is 
related to A by Aa = AaA. From the renormalized pa¬ 
rameters, the susceptibility is obtained accurately using 
renormalized perturbation theory (RPT) to second order 
in U/(ttA) M The result is 


the scope of the present paper. This precludes to give ex¬ 
act results for cb out of the symmetric case Ed = —U/2. 

In the symmetric case, V. Zlatic and B. Horvatic have 
shown that the exact results for the magnetic susceptibil¬ 
ity and quasiparticle weight z = A/A can be expressed 
as a convergent power seriesi^i We have evaluated this se¬ 
ries using a FORTRAN program in quadruple precision. 
From these results we have derived the renormalized pa¬ 
rameters and the Wilson ratio for comparison. 


X = (5AiB)Vd(0)R/2, 


( 8 ) 


D. Interpolative perturbative approximation 


where the Wilson ratio and the quasiparticle spectral 
density at the dot are 


Pd{oj) 


R = l + [/pd(0), 
A/tt 

(a; - + A 2 ’ 


(9) 
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Using this procedure we obtain for example for U = 
lOA and Ed = —U/2, a Wilson ratio R = 1.9939, while 
the exact result from Bethe ansatz (see Section IIII Cp 
is 1.9957 and NRG with the full density matrix gives 
2.024.^ 

The quasiparticle spectral density is modified by the 
renormalized interaction To calculate the width of 
the resulting renormalized spectral density we use ordi¬ 
nary perturbation theory (PT) ^^i^^ to second order in 
t//( 7 rA), taking Ed as the effective dot energy. Since 
even for U —>■ 00 , t//( 7 rA) is of the order of 1 (see Ta¬ 
ble HI, the second order results are accurate enough. To 
illustrate the consistency of the procedure, if one uses di¬ 
rectly PT for U = 3A in the symmetric case Ed = —U/2, 
one obtains a half width at half maximum Ap of the low- 
energy peak in the spectral density of the dot, only about 
2% higher (Ap = 0.558A) than the result Ap = 0.545A 
that comes determining first the renormalized parame¬ 
ters by NRG and then using PT. As Tq, Ap is also of the 
order of the Kondo temperature Tk in units for which 
the Boltzmann constant fcs = 1. 


C. Bethe ansatz 

As a check to the results obtained with the methods 
described above, we have calculated the occupancy and 
magnetic susceptibility for [/ —>■ 00 and several values of 
a shifted dot energy E^ at T = 0 using integral analyt¬ 
ical expressions obtained with the Bethe ansatz in Ref. 
[71I and as a particular case of a more general model in 
Ref. [ 72 I. Some tricks to deal with singularities in the in¬ 
tegrals were used, which are explained in the appendix of 
Ref. [ 7 ^ . Unfortunately, to obtain an analytical exact ex¬ 
pression for the second derivative of the occupancy with 
magnetic field d'^Ud/dB^, is very difficult and is beyond 


For small values of t// A one may use the interpolative 
perturbative approximation (IPA) which is a modifica¬ 
tion of the PT approach to second order in U/A modified 
to reproduce exactly the atomic limit U/A —^ -|-ooi2^^— 
In addition, the on-site term is split as ^dn-da = 
So- Ef^rida + J2ai^d “ Ef^)nda, where the second term 
is included in the perturbation, and is determined 
to optimize the results. In particular, in Ref. |^, A))® 
was determined to satisfy the Friedel sum rule. Results 
for the conductance through a QD using the IPA2^ agree 
with more recent ones using the finite temperature den¬ 
sity matrix renormalization group methodi^ 

Here we use the spin dependent version in which A^® 
depends on spin and is determined self-consistently for 
each spin to satisfy the corresponding Friedel sum rule 
Eq. ([2|i^ Gomparison with exact results without mag¬ 
netic field shows that the IPA provides very accurate val¬ 
ues for the occupancy Ud as a function of Ed for U/A < 6 . 
The maximum deviation is less than 1 % for U/ A = 61 ^ 


IV. RESULTS 

In this section, we provide results of the coefficient 
Cb of the magnetic-field dependence of the conductance 
[see Eq. ®] and the energy scale Tq as a function 
on the energy level of the QD Ed which is easily con¬ 
trolled by a gate voltage in transport experiments. Since 
To is determined by the magnetic susceptibility x [see 
Eq. ([SI] which is not accessible in transport experi¬ 
ments, we also provide a relation between Tq and the 
half width at half maximum of the spectral density Ap. 
In experiments with high asymmetry in the coupling be¬ 
tween the dot and both leads and in the Kondo regime, 
the differential conductance as a function of bias voltage 
G(I4) = dl/dVb directly represents the Kondo peak in 
the spectral densityi^i^ For less asymmetry and enter¬ 
ing in the intermediate-valence regime, the half width at 
half maximum of the zero-bias peak in G times the elec¬ 
tric charge e (Ag) increases to about 1.6 Ap but remains 
of the same order of magnitude for [/ —>■ ooi^ 

We focus our study in three values of U/A: 3, which 
is the largest value for which previous NRG results were 
reported^, U —>■ 00 which corresponds to the molecu¬ 
lar QDs, and an intermediate case U/A = 8 which is a 





reasonable value for semiconducting QDs in the Kondo 
regime. 


A. 17/A = 3 

In Fig. [T]we represent the coefficient of the magnetic- 
field dependence of the conductance [see Eq. dSD] as a 
function of the on-site energy for U = 3A and three dif¬ 
ferent methods. Because of a particular electron-hole 
symmetry (see Section II.C.l of Ref. HI), csiEd) = 
cb{—U — Ed), or in other words, for —U < Ed < —U/2, 
cs is the specular image around Ed = —U/2 of the re¬ 
sults represented in Fig. [TJ Our DMRG (NRG -I- RPT) 
results were obtained in intervals of A/4 (A/2) between 
the symmetric case Ed = —U/2 and Ed = 0. These 
results agree between them and also with the corrected 
NRG ones of Merker et al.i^ To compare with these re¬ 
sults, we have digitalized the results of given in Fig. 
1 of Ref. [S^ and used the relation 


cb \ Tq J 


X 

^sym 


2 
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where the superscript “sym” refers to the symmetric case 
Ed = —U/2. There is a small discrepancy for Ed = 
0, where the results of Merker et al. lie slightly above 
ours. This might be due to the fact that their results 
were calculated directly from the conductance, while ours 
use only thermodynamic quantities and are expected to 
be more accurate. In fact we have also calculated cb 
from the conductance derived from the spectral density 
calculated with the full-density-matrix NRG for a few 
points, but the results showed some deviations from the 
results obtained with DMRG and NRG-I-RPT using Eq. 

We also show in Fig. [T]the result of Cb within DMRG 
including only the first term of cb [neglecting the sec¬ 
ond negative term proportional to d^Ud/dB^ in Eq. ®] 
and com par ed with the corresponding result in Fig. 6 
of Ref. [2^ There is again a good agreement with 
the results of Ref. indicating a coincidence (ex¬ 
cept for deviations imperceptible in the figure) between 
the corresponding results for the occupancy Ud in the 
range of Ed studied. Although not shown in the fig¬ 
ure, the first term of cb (which depends only on the oc¬ 
cupancy Ud) within NRG-I-RPT and IPA also coincides 
with the results shown, because all these methods pro¬ 
vide accurately values for Ud- Instead the full result re¬ 
quires the calculation of the magnetic susceptibility x and 
d^Ud/dB'^ for which the IPA fails at large U as explained 
below. 

The correction due to the second derivative of the oc¬ 
cupancy becomes significant in the intermediate-valence 
regime, in particular ai Ed = 0, where it leads to a 
change of sign in cb [for all for all parameters consid¬ 
ered in this paper, the first (second) term in Eq. (6) is 
positive (negative)]. The effect of this correction (or in 



FIG. 1: (Color online) cb vs Ed for U = 3A obtained by 
different methods. 


general the difference between different approximations 
for Cb) is not so evident in the curve vs Ed because 
the magnetic susceptibility decreases strongly as Ed en¬ 
ters in the intermediate-valence regime (particularly for 
large U) and then the values of |c^| are largely reduced 
with respect to cb [by a factor 5.5 according to Eq. (fTTl) 
for Ed = Q]. The different approaches lead to the exact 
result Cb = c'^ = 7r^/I6 ~ 0.617 in the symmetric case 
Ed = — C//2^“—iSL and for large enough U the values of 
c'g lie near zero for Ed = 0, even neglecting d'^Ud/dB'^, 
because of the factor (Tq^™/To)^. 

The values of cb obtained using IPA lie slightly below 
those of the other methods (in vs Ed the maximum 
deviation in near 0.02). We have verified that the oc¬ 
cupancy is very well reproduced by the IPA (with an 
underestimation of the order of 1 %) in agreement with a 
previous comparison with Bethe-ansatz resultsi^ There¬ 
fore, the first term (positive in the range of Ed shown) of 
Cb in Eq. ([5]) is well reproduced. However, the remaining 
negative term is overestimated due to an underestimation 
of the magnetic susceptibility y (overestimation of Tq) 
by 6% and an overestimation of d^Ud/dB^ that reaches 
of the order of 10% for Ed near —A. The accuracy of 
the IPA depends strongly on the perturbation parameter 
U/{TrA). For example in the symmetric case Ed = —U/2, 
the IPA result for y is below the Bethe ansatz results by 
15% for U = 4A and by only 1.4% for U = 2A. 

One disadvantage of cb with respect to Cg is that the 
energy scale used as a reference, Tq, depends on Ed [see 
Eq. ([5])] and one needs to know two numbers for each Ed 
{cb and Tq) to describe the magnetic-field dependence of 
the conductance. An advantage is that Tq is of the or¬ 
der of the Kondo temperature, which in turn is propor¬ 
tional to the width 2AG/e of the zero-bias peak in the 
conductance (which is experimentally accessible) and to 
the width 2Ap of the spectral density of the dot state. 
All these quantities are of the order of the quasiparticle 
level width A. In Table U we give values obtained from 
RPT-I-NRG of these quantities, except Aq. This width 










is addressed in Ref. [13 where it is shown that in the 
Kondo regime for strongly asymmetric leads (as usual in 
molecular QDs) Aq = Ap.— 
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TABLE I: Effective parameters, half width of the spectral 
density Ap and ratio Ap/To obtained from NRG+RPT for 
U = 3A and several values of Ed- 


Ed/A 

A/A 

Ed/A 

U/{nA) 

Ap/A 

Ap/To 

-1.5 

0.639 

0 

0.738 

0.835 

0.944 

-1 

0.671 

0.196 

0.732 

0.848 

0.885 

-0.5 

0.754 

0.421 

0.716 

0.883 

0.766 

0 

0.845 

0.700 

0.698 

0.930 

0.581 

The 

Bethe-ansatz result 

of A for 

Ed = 

-1.5A is 


0.6522, slightly larger than the RPT+NRG result 0.639 
tabulated. This is due to the fact that for this calculation 
we used a finite half band width D = lOA that affects the 
Kondo temperature. However the Wilson ratio 1.738 is 
almost the same as the exact one 1.741. The ratio Ap/Tg 
evolves from near 1 in the symmetric case Ed = —U /2 to 
near 1/2 for Ed = 0- From the data given in Table |T] and 
Fermi-liquid relations one can obtain the occupancy 



FIG. 2: (Color online) occupancy and scaled magnetic sus¬ 
ceptibility X* = xA/[100((;/is)^] for 1/ —^ oo as a function of 
the shifted dot energy for different techniques. 
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which coincides in general within 1 % with the corre¬ 
sponding result obtained directly with NRG. 


B. U ^oo 

Experimentally, the on-site energy of the dot level Ed 
is controlled by the gate voltage and determined by ca¬ 
pacitance effects of all applied gate voltage s (see 
for example the supplementary material of Ref. 0 ). The 
position of the Goulomb blockade edges and the related 
charge-transfer peaks in the spectral density of the dot 
level are determined by a shifted energy E^ which con¬ 
tains the effects of a renormalization due to the hybridiza¬ 
tion with the leadsj ^^i^^ This shift is significant when 
both, the half band width D and the Coulomb repulsion 
U are much larger than A, as in many realistic systems. 
In particular for U —>- 00 , the calculation by Haldane 
based on poor man’s scaling gives 

6 = E*d-Ed = ^In , (13) 

where a ^ 1. Experimentally, E'^ is accessible but not 
Ed- For example at temperatures above the Kondo tem¬ 
perature, a maximum in the equilibrium conductance 
takes place for E^ = 0, or for finite E^ Coulomb edges 
appear for bias voltages such that Vb = EE^je^ 


In order to compare the DMRG results for the oc¬ 
cupancy Ud and magnetic susceptibility y for fixed Ed 
with the corresponding analytical Bethe-ansatz (in which 
ZI —>■ 00 is assumed) as a function of E^, we have to de¬ 
termine the shift (5 dmrg- We obtain that both functions 
Ud and X shifted by the same i5dmrg = 1.53A coincide 
with the corresponding Bethe-ansatz results, as shown 
in Fig. [21 This result is consistent with Eq. (IT^ which 
gives S = 1.47A for D = lOOA and a = 1. The same hap¬ 
pens with the NRG-bRPT results using (5 nrg = l.OOA. 
For large negative values oi Ed/A the occupancy flattens 
near 1, and the susceptibility increases strongly due to 
the exponential dependence of the Kondo temperature 
Tk on E^ and the fact that y is proportional to 1 /Tk in 
the Kondo regime. 

Once the shifts 5 are determined, we can represent cb 
as a function of E^ using both methods. The results are 
shown in Fig. |3l For E^/A < —2, as the occupancy Ud > 
0.84 and the system is in the Kondo regime, cb > 0.5 
indicating a small to moderate deviation from the value 
(7r/4)^ Ri 0.617 of the symmetric case. In this region, the 
correction due to the second (negative) term of cb in Eq. 
m is below 1 % Instead, for E^/A > —2, cb decreases 
rapidly and changes sign for E)J/A r; 0, in the middle of 
the intermediate-valence region. 

In Table El we list the renormalized parameters and 
ratios of different quantities proportional to Tk obtained 
with NRG-bRPT. As it is expected in the Kondo regime 
—Ed ^ A and Ed + U A, the Kondo temperature 
depends exponentially on Ed for Ed/A < —2. However, 
in this regime, particularly for Ed/ A — —3, the ratios 
Ap/A and Ap/Tg are almost constant, near 0.7 and 0.85 
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FIG. 3: (Color online) cb vs E’^ for 17 —^ oo obtained by 
DMRG and NRG+RPT. 


TABLE II: Effective parameters, half width of the spectral 
density Ap and ratio Ap/To obtained from NRG+RPT for 
U ^ oo and several values of Ed- 


Ed/A 

A/A 

Ed/A 

U/{-kA) 

Ap/A 

Ap/To 

-6 

2.51 X 10”'^ 

0.0937 

1.009 

0.706 

0.892 

-5 

1.21 X 10"^ 

0.118 

1.014 

0.705 

0.886 

-4 

5.79 X 10"® 

0.160 

1.025 

0.703 

0.873 

-3 

0.0270 

0.243 

1.054 

0.699 

0.838 

-2 

0.115 

0.416 

1.136 

0.693 

0.740 

-1 

0.356 

0.766 

1.317 

0.716 

0.530 

0 

0.640 

1.338 

1.594 

0.793 

0.286 


respectively. Instead, entering the intermediate valence 
regime, Ap/A slightly increases and Ap/To strongly de¬ 
creases. 

C. U/A = 8 

In Fig. 0] we show the cb for U = 8A obtained by 
DMRG and NRG+RPT. While the results of both meth¬ 
ods coincide in the Kondo regime {Ed < — 2A in the 
figure) there are some deviations in the intermediate va¬ 
lence regime. This is due to the fact that for finite U the 
Haldane shift 6 depends on Ed-, being 0 in the symmetric 
case Ed = —U/2, positive for Ed > —U/2 and larger for 
larger D- Therefore, the DMRG results calculated with 
D 10 times larger, correspond to larger in general, 
but for both techniques E*i = Ed in the symmetric case. 

Table IIIII displays the renormalized parameters and 
ratios Ap/A, Ap/To obtained with NRG+RPT. As in 
the case 17 —> oo, these ratios are nearly constant in 
the Kondo regime (although the variation is more pro¬ 
nounced as before) —Ed ^ A and Ed + U » A, while 



FIG. 4: (Color online) cb vs Ed U = 8A obtained by DMRG 
and NRG+RPT. 


TABLE III: Effective parameters, half width of the spectral 
density Ap and ratio Ap/To obtained from NRG+RPT for 
U = 8A and several values of Ed- 


Ed/A 

A/A 

Ed/A 

U/{nA) 

Ap/A 

Ap/To 

-4 

0.120 

0 

0.985 

0.715 

0.902 

-3 

0.143 

0.101 

0.987 

0.718 

0.895 

-2 

0.235 

0.247 

1.004 

0.724 

0.845 

-1 

0.457 

0.510 

1.040 

0.759 

0.699 

-0.5 

0.609 

0.515 

1.060 

0.802 

0.573 

0 

0.746 

0.977 

1.081 

0.857 

0.430 


in the intermediate valence regime, Ap/A increases and 
Ap/To decreases markedly (although not so strongly as 
for U -T oo). The Bethe ansatz result of A/A = 0.1326 
for Ed/ = —4A is about 10 % larger than the NRG+RPT 
result. As for U = 3A this difference is due to the effect 
of the finite band width 2D of the conduction band on 
the Kondo temperature, and practically does not alter 
the Wilson ratio. 


D. Cb as a function of occupancy 

In Fig. Owe represent our main results of cs as func¬ 
tions of the occupancy and the corresponding result for 
U ^ oo obtained previously using slave bosons in the 
mean-field approximation (SBMFA)i^ We also include 
the results using the interpolative perturbative approxi¬ 
mation (see Section lill Dl) for [7 = 2. In spite of the very 
different values of U used, all curves look qualitatively 
similar. There is a slightly more pronounced decrease of 
Cb vs 1 — Ud for smaller values of U, cb changes sign 
at Ud ^ 0.64 for [7 = 2 and at Ud ^ 0.57 for [7 —>• oo. 
Guriously, the SBMFA gives values that lie in between 






























FIG. 5: (Color online) cs vs 1 — Ud for different values of U 
obtained by DMRG and NRG+RPT, and for t/ —>■ oo using 
slave bosons in the mean-field approximation. 

those for U = 3A and those for U = 8A, although it is 
intended for 17 —>■ oo. 


V. SUMMARY 

Using different techniques (mainly DMRG and 
NRG-I-RPT) we have calculated the coefficient of the 
magnetic-field dependence of the conductance [see Eq. 
®] for several values of U focusing on the strong¬ 


coupling limit U ^ Tk- As it is known, cb = 
in the symmetric case = —C//2, where the occu¬ 
pancy Ud = 1. As Ed increases or decreases, cb decreases 
and changes sign in the intermediate valence cases where 
Ed ^ cp 01 Ed + U ^ ep, where E^ is a renormalized 
energy that includes the Haldane shift (which is impor¬ 
tant for large U) and ej’ is the Fermi energy. In these 
cases (related by a special electron-hole transformation) 
the occupancy is near 0.6 or 1.4 respectively. 

We have also provided quantitative results for the ra¬ 
tio of the energy scale of the magnetic-field dependence 
of the conductance and the magnetic susceptibility Tq 
[see Eq. ®] with the half width at half maximum of the 
spectral density Ap. For devices in which the coupling to 
the left and right leads are very different, this quantity 
coincides with the half width at half maximum Aq of the 
zero-bias peak in G times the electric charge e and is ex¬ 
perimentally accessible i^^’^° For more symmetric devices 

t /Ap increases and some values are tabulated in Ref. 

When Cb is represented as a function of the occupancy 
Ud for different values of U, the curves cpind) look qual¬ 
itatively similar. 
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